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Abstract 


-  y  The  paper  develops  a  theory  of  tearing  for  circumferential  through-cracks  in  pipes. 
Resistance  to  tearing  is  postulated  to  be  a  consequence  of  invariance  in  shape  at  the  tip  of  a 
growing  crack. 


NOTATION 


Young's  modulus 
fracture  parameter 
applied  bending  moment 
flow  stress 
M/(noFhR2) 


pipe  mean  radius 
pipe  thickness 
axial  displacement 
(h^)[12(l-v2)]-i/2 

E5/(GpR) 


INTRODUCTION 


A  method  of  J-estimation  for  a  circumferentially  through-cracked  pipe  subjected  to  bending 
loads  was  presented  in  [1].  The  method  was  based  on  the  Dugdale  model  and  the  semi-membrane 
theory  of  cylindrical  shells.  The  results  of  that  analysis  can  be  used  in  conjunction  with  a  material 
resistance  curve  to  examine  questions  of  crack  growth  stability  by  well-known  procedures.  The 


4  y 


’S'  'V 


’ *  v  v  v  .  V 


r  v ,  >* 


□a 


-2- 


present  paper  is  concerned  with  the  use  of  a  crack  growth  hypothesis  in  place  of  a  resistance  curve 
to  serve  as  an  alternative  theoretical  foundation  for  stability  analyses.  The  idea  for  such  an 
hypothesis  has  appeared  repeatedly  in  the  literature  [2]  -  [7]  in  the  form  of  a  statement  that  some 
quantity  or  other  associated  with  the  crack  tip  remains  constant  as  the  crack  grows.  A  common 
feature  of  such  crack  growth  or  "tearing"  hypotheses  is  that  they  lead  to  a  nonlinear  first-order 
differential  equation  for  J  as  a  function  of  crack  extension.  In  other  words  the  resistance  curve 
itself  is  determined  theoretically  instead  of  experimentally.  The  determination  of  a  resistance  curve 
is  not  really  necessary  in  applications  of  the  idea  to  cracks  in  structures  as  the  present  paper  will 
show.  Such  an  approach  has  theoretical  appeal  partly  because  the  difficulties  with  geometry 
dependance  in  experimentally  determined  resistance  curves  are  avoided.  In  applications  the  effects 
of  geometry  dependance  are  an  essential  part  of  the  problem.  Of  perhaps  greater  theoretical  appeal 
is  the  prospect  that  a  growth  hypothesis  (provided  the  idea  proves  to  be  correct)  brings 
phenomenological  fracture  mechanics  one  step  closer  to  what  must  be  its  physical  foundations. 

The  tearing  hypothesis  used  in  the  present  paper  is  an  adaptation  to  the  pipe  problem  of  a 
similar  criterion  set  forth  in  [6]  for  the  case  of  cracks  in  flat  plates  according  to  a  Dugdale  model. 
The  analysis  is  based  on  results  from  [1]  and  leads  to  a  first  order  nonlinear  differential  equation 
governing  the  effects  of  crack  growth.  The  applied  bending  moment  as  a  function  of  the  crack 
advance  angle  and  the  increased  flexibility  of  the  pipe  due  to  the  crack  follow  directly  from  the 
solution  to  the  equation. 


ANALYSIS 

Figure  1(a)  shows  a  sketch  of  the  upper  half  of  the  crack  profile  at  the  instant  of  crack 
growth  initiation.  The  shaded  region  indicates  plastically  deformed  material  according  to  the 
Dugdale  model.  The  crack  tip  lies  at  0  *  a„ ,  and  the  Dugdale  zone  ends  at  9  =  p0 .  The  value  of 
u  at  0  =  a,,  is  half  the  critical  crack  opening  displacement  5C  related  to  the  critical  value  of  the 
fracture  parameter  Jc  and  the  flow  stress  by  the  familiar  relation 
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The  picture  after  the  crack  tip  has  advanced  to  0  -  a  is  shown  in  Fig.  1  (b).  An  hypothesis  of  the 
model  is  that  the  crack  advances  into  the  plastic  region  leaving  behind  a  strip  of  "dead"  material 
(shown  unshaded).  The  crack  face  is  thought  of  as  lying  along  the  lower  edge  of  the  dead 
material.  If  the  upper  curve  is  given  by  the  function  u(0,  a)  then  the  width  of  the  dead  strip  is  by 
definition  u(0, 6) .  The  width  of  the  crack  opening  is  given  by 

A  =  2E(0,a)- 211(0,0)  ,  a0  <  0  <  a  (2) 

The  usual  terminology  isn't  quite  appropriate  but  the  crack  tip  opening  angle  co  will  be  defined  to 


Rfi>  *  -  lim 

0— kx  <*0 


d5  Jd u 


where  5  =  2u(a,  a) .  With  use  of  J  =  Op5 ,  and  put  in  dimensionless  form,  Eq.  (3)  reads 


where 


E  dJ  .  c  E  r 


F--2 


Equation  (4)  might  be  called  the  "tearing"equatk». 

The  crack  growth  hypothesis  set  forth  in  [6]  is  that  the  crack  opening  shape  in  the 
neighborhood  of  the  crack  tip  remains  invariant  as  the  crack  grows.  The  obvious  adaptation  of  that 
hypothesis  in  the  present  case  is  that  the  angle  (0  remain  constant.  In  essence  the  hypotheses 
adopted  here  are  the  same  as  those  in  [6]  but  there  is  some  difference  in  detail  In  [6]  the  analysis 
is  made  according  to  plane  stress  theory  for  cracks  in  flat  plates.  Here  the  semi-membrane  theory 
for  cylindrical  shells  applies.  Results  from  semi-membrane  theory  are  certain  limits  of  results  from 
a  more  complete  shell  theory.  What  that  means,  for  instance,  is  that  the  angle  <5  given  by  (3)  is  a 
parameter  associated  with  the  shape  of  the  crack  opening  near  the  tip,  not  an  actual  angle. 
However,  it  is  a  somewhat  remarkable  fact  that  quantities  important  to  fracture  mechanics  turn  out 
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to  be  given  accurately  by  semi-membrane  theory  provided  the  crack  is  long  enough.  See  the 
discussion  after  Eq.  (62)  in  [1]. 

What  remains  to  be  done  is  to  relate  the  angle  dj  to  material  properties  and  the  radius  and 
thickness  of  the  pipe.  In  order  for  the  theory  to  make  sense  <5  cannot  in  any  way  depend  upon  the 
crack  length.  A  material  resistance  curve  is  generally  assumed  to  have  been  obtained  as  though 
small  scale  yielding  conditions  were  adhered  to.  A  crack  growth  hypothesis  such  as  the  one 
adopted  here  is  tantamount  to  assuming  that  the  whole  resistance  curve  is  determined  by  conditions 
at  the  initiation  of  crack  growth.  There  is,  however,  no  reason  to  suppose  that  55  cannot  depend 
upon  the  curvature  and  thickness  of  the  pipe.  Obviously  the  effect  of  curvature  on  55  must  be 
established  on  theoretical  grounds  because  curvature  is  not  a  material  property. 

The  required  determination  of  Q  is  accomplished  here  by  analysis  of  certain  hypothetical 
experiments.  Small  scale  yielding  conditions  must  exist  in  the  limiting  case  aP  ->  ®o  .  In  such  a 


case  the  extent  of  the  Dugdale  zone  is  small  and  results  from  [1],  repeated  here 

e6- V20-a)2 


eF  *  2V2  (p-a) 


imply 


in  which  case  the  tearing  equation  reads 


F-2 


5/4/  EJ 

^EOpR 


E  dJ  ,5/4  /  EJ  E 

tt  t:+2  t — rr-  =  t=t  Q 


OpR 


EOpR 


The  coefficients  in  the  equation  do  not  depend  on  a  therefore  the  solution  depends  only  on 
a-cxo-  Initially  J  -  Jc  and  initially  the  first  term  in  (9)  is  the  Paris  tearing  modulus  T0  obtained 
from  the  initial  slope  of  a  resistance  curve  for  small  scale  yielding.  The  angle  co  is  thus 
determined  to  be  given  by 


A  somewhat  different,  but  related,  calculation  was  made  based  on  the  analysis  of  a  semi-infinite 
circumferential  crack  in  a  shallow  cylindrical  shell.  The  mathematical  formulation  of  the  problem 
in  terms  of  a  semi-membrane  theory  turns  out  to  be  well-set  and  the  problem  can  be  solved.  The 
details  will  not  be  reproduced  here  but  the  resulting  determination  of  m  is  exactly  the  same  as 
(10).  Since  the  equations  of  the  model  are  indifferent  to  the  magnitude  of  the  parameters  involved 
the  determination  (10)  is  expected  to  hold  in  all  cases.  In  the  general  case  the  tearing  equation  has 
the  form 

in  which  F  is  a  complicated  implicit  function  of  J  and  a. 

The  following  expression  for  the  dimensionless  axial  displacement  on  the  cracked  section 
was  derived  in  [1] 

V2eu(0,a)  «Ccos0 -A  +  (3cos0--j0sinG)a  ,  |0|<a  (12) 

where  A  and  C  are  functions  of  a ,  {) ,  y,  and  c  (results  are  reproduced  in  the  Appendix).  The 

load  parameter  o  (proportional  to  the  applied  bending  moment  on  the  pipe)  is  also  a  function  of 

the  three  angles.  The  angle  y  marics  the  beginning  of  a  compressive  Dugdale  zone  on  die  side  of 

the  pipe  away  from  the  crack  (back  yielding).  In  the  absence  of  such  a  zone  y  =  x .  From  (12) 

and  definitions  the  tearing  equation  (1 1)  can  be  put  in  the  following  form 

C'cos  a-A'  +  (3cosa-^a  sina)c' «  to 

(13) 

where  co«(Ee/V2  Op)5 

and  where  a  prime  denotes  total  differentiation  with  respect  to  a  .  In  terms  of  A,  C,  a,  and  a 
the  fracture  parameter  J  is  given  by 

J*(V2  0pR/Ee)[Ccosa-A  +  (3cosa-^a  sin  a) a)  (14) 

Since  J  is  not  a  convenient  dependent  variable  f)  will  be  used  instead.  For  values  of  a 
less  than  that  for  which  back  yielding  appears  A,  C,  and  o  are  known  explicitly  in  terms  of  a  and 
B  alone.  Beyond  that  range  y  is  known  implicitly  as  a  function  of  a  and  3  and  a  root-finding 


routine  must  be  included  as  part  of  the  numerical  integration  procedure.  In  either  range  the 
differential  equation  for  p  has  the  form 

P'  =  f(P,a)  (15) 

The  initial  value  of  {3  is  determined  from  (14)  by  the  condition  J  =  Jc  at  a  =  Oq  ,  and  by  a  root¬ 
finding  process.  For  very  tough  materials  back  yielding  can  occur  before  crack  growth  initiation, 
in  which  case  initial  values  of  P  and  y  must  be  obtained  by  root  finding. 

The  procedure  for  setting  up  (15)  will  be  described  in  some  detail  in  the  case  y  <  it .  From 
results  in  [1]  there  are  two  expressions  for  o  involving  the  three  angles  a,  (3,  and  y.  These 
have  the  form 

N1(a,p,y)-oDl(P,Y)=0  (16) 

N2(a,p,Y)-oD2(p,Y)  =  0  (17) 


Eliminate  a  between  these  two  equations  to  obtain 


Ni  °2  "  N2  =  0 


(18) 


For  a  and  P  given  this  is  an  implicit  equation  for  the  determination  of  Y-  It  must  be  solved  by 
some  root-finding  method.  Then  o  is  given  by  Nj/Di  or  by  N2/D2  (which  are  equal).  The 
total  derivative  of  (16)  with  respect  to  a  yields 


-Dj  o’  =  - 


A  similar  equation  follows  from  (17)  containing  subscripts  2  in  place  of  1.  These  are  two 
simultaneous  equations  to  solve  for  y  '  and  o'  in  terms  of  P'.  Let  the  result  be  expressed  in  the 
form 


r=r1+r2p’  (20) 

o’-VIjp’  (21) 


Equation  (13)  can  now  be  written  in  die  form 

g2p,+  gi  *© 


(22) 


(23) 
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and  where 

q  =  3  cos  a  -  -i  a  sin  a  (25) 

The  f  in  (15)  can  now  be  expressed  as 

f«(a>-giyg2  (26) 

The  various  functions  required  to  calculate  gj  and  g2  are  given  in  the  Appendix. 

So  long  as  the  dimensionless  axial  stress  at  6  =  x  is  greater  than  minus  one  back  yielding 
does  not  occur.  This  condition  can  be  written 

o+  ^  V2  (P  -  a  -  a  sin  pysin  x  <  1 

(27) 

where  x  ■  (x  -  £)/V2 

Let  aB  denote  the  value  of  a  for  which  the  condition  is  first  violated.  For  a  <  aB  the  above 

formulas  simplify  because  y  =  x .  Equation  (16)  determines  a,  and  o'  is  obtained  from  (19)  by 
dropping  the  term  containing  y ' .  The  terms  with  Tj  and  T2  drop  out  of  (23)  and  (24)  and  the 

formulas  in  the  Appendix  hold  with  y  *  x .  The  numerical  integration  of  (15)  proceeds  in  two 
stages,  one  for  a  <  aB  and  another  for  a  >  aB .  Two  difficulties  arise.  One  is  simply  a  matter  of 
matching  the  two  stages.  In  any  step-by-step  process  a  is  not  likely  to  equal  aB  exactly  at  the  end 

of  a  step.  Some  sort  of  interpolation  scheme  is  required  to  ensure  accuracy  in  the  starting  values  of 
the  second  stage.  A  more  troublesome  difficulty  arises  from  the  fact  that  y  ■  x  is  always  a  root, 
but  not  the  desired  root,  of  (18).  the  spurious  root  can  be  eliminated  by  division  of  both  N2  and 
D2  by  x  -  y  and  using  a  polynomial  approximation  to  sin(z)/z  to  be  found  in  [9].  However, 
finding  the  root  when  x  -  y  is  small  is  still  difficult  Fortunately  it  turns  out  that  inaccuracies  in 
small  values  of  x  -  y  matter  very  little  later  on  in  the  calculations. 
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For  any  given  value  of  a  there  is  a  value  of  the  load  parameter  o  sufficiently  large  to  cause 

plastic  collapse.  This  limit  load  is  given  by 

0^=  J(cos|_|  sina)  (28) 

Some  example  results  are  shown  in  the  figures.  These  come  from  the  following  round- 
number  choices  for  the  parameters:  R/h  =  10,  oF  =  100  ksi,  E  =  30x1 06  psi,  Jc  =  1000  lb/in., 
R  =  15  in.,  T0  =  100  or  25.  Figure  2  shows  load  curves  for  initial  crack  angles  of  .3,  .5,  and  .7 
radians;  the  dashed  curve  is  for  T0  =  25.  The  points  where  a  -  aB  are  indicated  by  tick  marks. 
The  upper  curve  represents  the  limit  load  according  to  (28).  In  the  cases  where  T0  =  100  the  load 

reaches  a  maximum  slightly  below  the  limit  load  after  a  small  amount  of  tearing.  In  the  case 
T0  =  25  the  load  reaches  a  maximum  considerably  below  the  limit  load  after  a  larger  amount  of 

tearing.  Theoretically  the  load  curves  should  approach  the  limit  curve  asymptotically.  Numerically 
they  appear  to  become  tangent  to  the  limit  curve.  Typically  the  computer  indicates  an  overflow 
condition  and  the  program  may  as  well  be  stopped.  Figure  3  shows  computed  resistance  curves  in 
three  cases.  In  these  cases  (and  in  many  computed  but  not  shown)  the  initial  slope  corresponds 
very  closely  to  T0.  There  is  very  little  dependance  on  in  the  case  T0  =  100. 

CONCLUDING  REMARKS 

The  results  from  the  theory  appear  to  be  qualitatively  correct  but  no  extensive  comparisons 
with  the  results  of  experiments  have  yet  been  made.  A  tearing  theory  such  as  that  presented  here 
substitutes  a  crack  growth  hypothesis  for  the  hypothesis  that  there  exists  a  "material"  resistance 
curve.  The  theory  calls  for  the  experimental  determination  of  three  material  constants  Jc  ,  and  T 
instead  of  Jc  ,  and  a  resistance  curve.  Such  simplicity  is  appealing.  The  very  need  for  a  material 
resistance  curve  (always  a  rather  troublesome  object)  is  bypassed.  Should  the  theory  turn  out  to 
have  validity  a  search  for  a  more  realistic  model  accounting  for  strain  hardening  would  be 
indicated. 

As  previously  noted  the  numerical  implementation  of  the  theory  given  here  runs  into 
difficulties  when  the  load  parameter  gets  close  to  the  plastic  collapse  value.  Examination  of  the 
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limiting  case  P  -  y  small  or  even  P  -  y  -  0  (fully  yielded  cracked  section)  has  not  yet  been  carried 
out  analytically.  Because  shell  theory  is  involved  results  in  the  limiting  case  should  be  similar  to 
but  not  identical  to  results  from  theories  such  as  those  in  [7]  or  [8]  in  which  the  cracked  section  is 
assumed  to  be  fully  yielded  from  the  outset.  The  crack  growth  hypothesis  in  [7]  bears  an 
interesting  similarity  to  the  one  adopted  in  the  present  paper. 
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APPENDIX 


Q  =  a  sin  y-  (n  -  y) 

P  sin  x  =  Q  cos  x  -  a  sin  P+  p  -  a 


where  x  =  (y  -  P)/V2 

A  =  ocos p+  ^ (P  -  a)2-  V2  (P cos  x  +  Q  sin  x) 

C  sin  p  =  -  ^  (7  sin  P+ pcosP)a+ P  -  a  +  sin  (P~a) 


^  =  -(P-a)  +  V2  cot  x 


=  -  (P  -  a  -  a  sin  P)cot2  x  -  V2  (1  -  o  cos  P )cot  x  -  [a  sin  y  -  (n  -  y)]cot  x/sin : 
dp 


^  -  -  V2  (1  +  a  cos  y)/sin  x  +  {[osin  y-  (jc  -  y)]cos  x  +  p-a-asin  p}/sin2  x 


^  =  cos  P  -  V2  sin  y/sin  x  +  V2  sin  P  cot  x 


^  -  [1  +  cos  (p  -  a)]/sin  P 
da 


^  =  [^  (P -  sin  P  cos  p)  0+  sin  a  +  sin  P  -  (P  -  a)cos  P]/sin2 P 


dC 


3? 


=  0 


^  1  (7  sin  P  +  p  cos  P)/sin  P 


Nj  =  [sin  P  -sin  a+  (P  -  a)cosP]sin  x+  V2  (p  -  a)sin  P  cos  a  -  V2  (tt  - y)sin  P 
=  -  (cos  a  +  cos  P)sin  x  -  V2  sin  p  cos  x 


=  2  cos  P  sin  x  +  -i  V2  [sin  P  +  sin  a  +  (P  -  a) cos  P]cos  x-  V2  (7t  -  y)cos  P 
dp  ^ 


aNj 


-=^  =  ^  V2[sin  P  -  sin  a  +  (P  -  a)cos  p]cos  x-  (P  -  a)  sin  p  sin  x  +  V2  sin  P 
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D^  =  I(p+3sinPcosP)sinx+  V2  sin2  P  cos  x-  V2  sin  y  sin  P 

B.o 

da 

^.  =  2  cos2p  sin  x- jV2(p-5  sin P  cosp)cos  x- V2  sin  y  cos  p 
^-  =  ^V2(p+ 3  sin  P  cosp)cos  x  -sin2p  sin  x  -  V2  cos  y  sin  P 

N2  =  [sin  y-(jc  -y)cosy]sin  x  +  V2(ji  - y);in  ycosx-V2(P  -a)sin y 

^2  -  V2  sin  y 
9a 

=  -jV2 [sin y - (jc - y) cos  y]cos x  +  (« - y) sin y  sin x - V2  sin y 

^2.  =  2  cos  y  sin  y  -  A  V2  [sin  y  -  (n  -  y)  cos  y]cos  x  -  V2  (P  -  a)  cos  y 
9y  2 

-  2  (jt  - y-  3  sin  y  cos  y)sin  x+  V2  sin2y  cos x-  V2  sin  y  sin  P 

H‘° 

a-2.V2(Jt-y-3sinycos y )cos x  +  sin2 y  sin x - V2  sin y  cos  P 

2^  =  -2  cos2  y  sin  x  +  j  (it  -  y  +  5  sin  y  cos  y)cos  x  -  V2  cos  y  sin  p 
dy  4 
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